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Abstract: We compute the phase diagram of the quark- meson model at finite temperature, finite 
baryon chemical potential [lb = 3/i and constant external magnetic field B, using the functional 
renormalization group. Our results show that the critical temperature increases as a function of B 
at fi = 0, but for values [i larger than about 210 — 225 MeV, the opposite behavior is realized. As 
the magnetic field increases, the critical point (T*,/x*) moves from large fi, small T towards small 
/i, larger T in the [i-T phase diagram. 
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1 Introduction 

The phase structure of QCD as a function of temperature T and baryon chemical \xb has received 
a lot of attention. This interest has been spurred by the ongoing heavy-ion collision experiments at 
RHIC (Brookhaven) and the LHC (CERN) and the search for the location of the critical endpoint 
where the curve of first-order chiral transitions terminates. Moreover, it turns out that the phase 
diagram has a much richer structure than was anticipated. In particular, there are a number of color 
superconducting phases at large baryon chemical potential and low temperature (see for example 
the reviews [1, 2] for a thorough discussion). More recently, there has been quite some interest in the 
behavior of strongly interacting matter in external electromagnetic fields, motivated by astrophysical 
phenomena as well as heavy-ion collisions. Large magnetic fields exist inside ordinary neutron stars 
as well as magnetars [3]. Similarly, strong, but short-lived magnetic fields, up to B ~ 10 19 Gauss 
or \qB\ ~ 6m^, may be generated during noncentral heavy-ion collisions at RHIC and LHC [4-6]. 
In the context of strongly interacting matter in an external magnetic field, at least three important 
questions have arisen: Does the order of the transition change with the strength of the magnetic 
field B?; how does the critical temperature T c depend on the strength of the magnetic field B?; 
and is there a splitting between the chiral transition and the deconfinement transition beyond some 
finite value of the I?-field? We will address the first two of these questions in the present work. 

In Ref. [7] chiral perturbation theory at leading order was used to study the chiral transition as a 
function of magnetic field. Comparing the pressure of a hot pion gas with that of a noninter acting 
quark-gluon plasma with a subtracted vacuum energy term due to a nonzero gluon condensate, they 
found a first-order transition for weak magnetic fields. The line of first-order transitions ends at a 
critical point (y/\qB\,T) = (600,104) MeV, where | q\ is the pion electric charge. For larger values 
of \qB\, the transition is a crossover. 

The authors of Ref. [8] use the Polyakov loop extended quark-meson (QM) model to study the 
phase transition. The bosons are treated at tree level while the fermions are treated at the one-loop 
level, except that the vacuum fluctuations at zero magnetic field are neglected. The authors report 



that they find a first-order transition at the physical point if they keep the remaining S-dependent 
vacuum fluctuations and a crossover if they are ignored. The mean-field calculations of Ref. [9] as 
well as the functional renormalization group calculation of Ref. [10] using the QM model suggest 
that the transition remains a crossover. The conflicting result between Ref. [8] and Refs. [9, 10] 
could be due to a different treatment of the vacuum fluctuations. In Refs. [11, 12] the effects of the 
vacuum fluctuations on the chiral transition were studied for B = 0, and it was shown that both the 
critical temperature T c and the order of transition depend on whether they are included or not. 

Most model calculations conclude that the critical temperature is an increasing function of the 
magnetic field B. This includes both mean-field type calculations [8, 13-20] involving (Polyakov- 
loop extended) Nambu-Jona-Lasinio (NJL) and quark-mesons models as well as chiral perturbation 
theory, and beyond mean field using functional renormalization group techniques [10]. These results 
are in line with the lattice calculations of Refs. [22, 23], however, these were carried out for bare 
quark masses that correspond to very large pion masses in the 200-480 MeV range. The very recent 
simulations of Bali et al [24, 25] for a pion mass of 140 MeV as well the MIT-bag calculations of 
Ref. [26] suggest that the critical temperature decreases with the strength of the magnetic field. 

In order to address the issue of whether there is a splitting between the chiral and deconfinement 
transitions, one must go beyond the NJL or QM models since these have nothing to say about 
the latter. This can be done by introducing an effective phenomenological potential for a constant 
temporal gluon field and expressing it in terms of the thermal expectation value of the trace of 
the Polyakov loop [27] . The conclusion of Ref. [8] is that the deconfinement transition temperature 
decreases with the strength of the magnetic field, while the chiral transition temperature increases. 
Similar conclusions were reached in [10], where a functional renormalization-group approach was 
used together with the P(QM) model to include mesonic fluctuations in the calculations. 

In the present paper, we continue to explore the properties of QCD in an external magnetic field 
using the quark-meson model. The work is a continuation of Ref. [9] in which the phase diagram 
in the \x~T plane was mapped out for strong fields, \qB\ ~ 5m^. The calculations were carried out 
in the mean-field approximation and the large- iV c limit where the bosonic modes are treated at tree 
level. Moreover, the influence of the fermionic vacuum fluctuations on the transition was studied. 
In the chiral limit, it was found that the transition is first order in the entire \i — T plane if vacuum 
fluctuations are not included and second order if they are included. Here we are using the functional 
renormalization group (FRG) [28] to investigate the behavior of strongly interacting matter in a 
magnetic background. The functional renormalization group is a non-perturbative, self-consistent 
resummation, with a wide range of applicability. See for example [30] for a recent review. 

The paper is organized as follows. In Sec. 2, we briefly discuss the quark-meson model and the 
functional renormalization group. In Sec. 3, we present and discuss our numerical results while we 
summarize in Sec. 4. The renormalization group equation for the effective potential is derived in 
Appendix A. 



2 Quark-meson model and the functional renormalization group 



The 0(4)-invariant Euclidean Lagrangian for the quark- meson model is 
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where a is the sigma field, 7r denotes the neutral and charged pions. r are the Pauli matrices, 
I 1 = + A*d) is the quark chemical potential, in terms of \i u and the chemical potential for the 
u and d quarks, respectively. The baryon chemical potential is given by \xb = 3/i. We set [i u = fid 
so that we are working at zero isospin chemical potential, /ij = ^(fi u — fid) = 0. The Euclidean 7 
matrices are given by 7^ = 27^, 74 = 7^, and 75 = — 7^, where the index M denotes Minkowski 
space. The fermion field is an isospin doublet 



(2.2) 



If h = 0, Eq. (2.1) is invariant under 0(4). If h ^ 0, chiral symmetry is explicitly broken, otherwise 
it is spontaneously broken in the vacuum. Either way, the symmetry is reduced to 0(3). Note that 
this requires m? < which is assumed in the remainder of the paper. 

Chiral symmetry is broken in the vacuum by a nonzero expectation value <f> for the sigma field 
and we make the replacement 



a 



where a is a quantum fluctuating field. The tree-level potential then becomes 

U A = + ^ - H ■ 



(2.3) 



(2.4) 



Note that we have introduced a subscript A on U, m 2 , and A, where A is the ultraviolet cutoff of 
the theory (see section 3). This is a reminder that these are unrenormalized quantities 1 . 

We will follow Wetterich's implementation of the renormalization group ideas based on the effective 
average action [(f)] [28] . This action is a functional of a set of background fields that are denoted 
by (f>. Tk[</>] satisfies an (integro-differential) flow equation in the variable k, to be specified below. 
The subscript k indicates that all the modes p between the ultraviolet cutoff A of the theory and 
k have been integrated out. When k = A no modes have been integrated out and Ta equals the 
classical action S. On the other hand, when k = 0, all the momentum modes have been integrated 
out and Tq equals the full quantum effective action. The flow equation then describes the "flow" in 
the space of effective actions as a function of k. 

In order to implement the renormalization group ideas, one introduces a regulator function Rf.(p). 
The function Rk{p) is large for p < k and small for p > k whenever < k < A, and R\(p) = 00. 



1 The symmetry breaking term is equivalent to an external field that does not flow and therefore h — /ia- 



These properties ensure that the modes below k are heavy and decouple, and only the modes 
between k and the UV cutoff A are light and integrated out. The choice of regulator function has 
been discussed in detail in the literature and some choices are better than others due both to their 
analytical and stability properties .We briefly discuss our choice of regulator in Appendix A. 

The flow equation for the effective action cannot be solved exactly so one must make tractable 
and yet physically sound approximations. The first approximation in a derivative expansion is the 
local-potential approximation and in this case the flow equation for reduces to a flow equation 
for an effective potential Uk(<p)- In the case of a constant magnetic field, the differential equation 
for Uk reads [10] 

]A ( i i 

d k U k = { [1 + 2n B (ui 1>k )] + [1 + 2n B (u 2 ,k)] 




where we have defined U\ ik = \/k 2 + U' k , W2,fc — \/k 2 + TP + 2U' k f p , u q> k = \/k 2 + 2g 2 p , p\(q, m, s) = 
(2m + 1 - s)\qB\ , n B (x) = l/(eP* - 1) , n$,(x) = l/(e^ x± ^ + 1), and p = \4> 2 . We briefly discuss 
the derivation of Eq. (2.5) in Appendix A. If we neglect the bosonic fluctuations by ignoring the 
first two lines of Eq. (2.5), we can solve the equation for U k analytically. This yields the standard 
mean-field result for the one-loop thermodynamic potential [9]. 

At zero temperature, the Bose distribution function vanishes and the Fermi distribution function 
becomes a step function. The flow equation (2.5) then reduces to 



dkU k = 




Furthermore, if we set pi = 0, the step function in (2.6) vanishes and we obtain the flow equation in 
the vacuum. 

3 Numerical results 

Solving the flow equation can in practice be done in two ways. One method is to make a polynomial 
expansion around the fc-dependent minimum and truncate the series, which leads to a set of coupled 
equations for the coefficients of the terms in this series. The other possibility is to solve the flow 
equation numerically for the potential Uk(4>) discretizing the field value cj>. We will do the latter, 
but discretizing instead in the quantity p = 4- in the interval p £ [0;8000] (MeV) 2 , using n = 200 




Figure 1. Potential Uk(4>) for different values of k. Note that the curves do not evolve monotonically in k. 
The inset shows the minimum (j>o = as a function of t = log ^. 



grid points 2 . The probed interval corresponds to ^ € [0; 126.5] MeV, which comfortably includes 
the minimum of the potential. 

Earlier we pointed out that h = h\ and it does not appear in the flow equation (2.5). However, 
one can show that h = f^vn^ and so a nonzero h is introduced via the boundary condition (2.4) of 
the effective potential. One must therefore tune the parameters m\ and Aa separately for h = and 
h 7^ to reproduce pion decay constant f n and the pion mass m n in the vacuum correctly. After 
tuning the parameters, the tree level potential U\{4>) is given by Eq. (2.4), with m\ = throughout, 
and Aa = 1-72 for 3 h = f-wirn^^ and Aa = 2.125 for 4 h = 0. This ensures in both cases that at 
T = \i = B = 0, the minimum of the potential is at 0n — fn — 93 MeV, and that the pion mass is 
correctly m n = 140 MeV for h > 0, h = 0, respectively. Also the a mass is in the allowed range of 
the broad resonance m a = 400 — 800 MeV, and we specifically take m a = 450, 480 MeV, again for 
h > and h = 0. We use an ultraviolet cutoff of A = 500 MeV as in Ref. [31]. We checked that 
changing the cutoff to 800 MeV results in a correction of approximately 3% to T c at \x = 0, B = 0, 
a precision we expect to persist throughout. The bare coupling of course changes significantly, but 
the physical result does not. 

We ignore the running of the Yukawa coupling and set g\. = 3.2258 for all k. This yields a quark 
mass of m q = gcfto = gf n = 300 MeV. The ^-dependent masses m 2 k and rr? a k are related to the 
/c-dependent minumum ^> 0j fc = fw,k of the effective potential as follows 

(3.1) 

<t> = fTT,k 

2 We checked that the results were insensitive to using half as many point, n = 100. 
3 We will call this the "h > 0" case or physical point. 
4 h — is the chiral limit. 
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Figure 2. Real part of the effective potential Uk(4>) for fi = T = B = 0. At tree-level Ua{4>) (black) and 
when fully renormalized Uk=o{<p) (green). The blue and red curves are the effective potential where we have 
ignored the effects of the bosonic quantum fluctuations and fermionic vacuum fluctuations, respectively. 

The physical masses m 2 and m 2 are given by evaluating Eq. (3.1) at k = 0. Finally, when we 
integrate the flow equation, we use the dimensionless variable t = log Thus k = A corresponds to 
t = and k = corresponds to t = — oo. 

In Fig. 1, we show how the potential Uk(4>) develops for different values of k. The black curve 
is the tree-level potential, k = A, and the green the full quantum effective potential Uo(tfi). It is 
interesting to note that the potential at intermediate stages, here shown as red, blue, and orange 
curves, do not evolve monotonically. This shows that the bosonic and fermionic terms in the flow 
equation dominate in different regions of integration. It is important to point out that the exact 
effective potential is both real and convex. It is an artefact of our truncation that the potential 
becomes nonconvex and that it develops an imaginary part [32]. 

The inset shows the minimum (j)Q = f n ^ as a function of the dimensionless variable t. The 
minimum converges to f n around t = —3, but other features of the potential, such as the curvature 
in the potential minimum, continue to develop somewhat later. By t = —6 the potential has settled. 

In Fig. 2, we show various instances of the quantum effective potential Uk(4>) for fj, = T = B = 0, 
h > 0, i.e. in the vacuum. The black line is the tree level potential Uk=A(<j>) and the green line the 
fully renormalized effective potential Uk=o(4>)- The red/blue lines are also fully renormalized, but 
ignoring the fermion/boson contribution to the fluctuations. We see that the two pull in opposite 
directions, enhancing/ weakening symmetry breaking. The difference comes basically down to a 
different sign of the relevant terms in the flow equation (2.6). In some cases, the fermionic quantum 
fluctuations can destabilize the effective potential altogether. Notice also that the green line is not 
simply the average of the red and the blue, as the flow equation performs non-linear resummations. 

Then in Fig. 3, we show the effective potential for different values of the magnetic field B measured 
in units of the ultraviolet cutoff A 2 , at T = fi = 0, and here shown for h > 0. We notice that the 
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Figure 3. Full quantum effective potential Uo{4>) for different values of the magnetic field. The arrow 
indicates the direction of increasing B. Inset shows the minimum (j)Q as a function of the magnetic field B 
with error bars. 

effective potential becomes deeper and the minimum moves to the right with increasing B. There 
is some roughness of the effective potential for p to the left of the minimum. The reason is that the 
potential is complex to the left of the minimum and numerical integration therefore becomes more 
noisy. We checked that this has no influence on the location of the minimum of the potential. The 
inset shows the minimum (f>Q as a function of the strength of the magnetic field. The minimum is 
within discretization errors an increasing, but non- linear, function of B. The error bars are meant 
to represent the resolution of the discretization in (f>. We found similar behaviour for h = 0. This 
result is in agreement with magnetic catalysis of dynamical symmetry breaking (MCDS) which is the 
effect that chiral symmetry is broken dynamically for any nonzero magnetic field when it is intact for 
B = 0, and more generally that a nonzero magnetic enhances symmetry breaking. This effect is now 
well established in fermionic systems in an external magnetic field, see for example Refs. [33-39]. 

In the left-hand frame of Fig. 4, we show the minimum of the effective potential for B = \i = 0, 
h > as a function of temperature T, allowing for the determination of the critical temperature at 
this value of the chemical potential. The order parameter goes to zero in a continuous manner which 
shows that the transition is second order or a cross-over. The finite coarseness of the discretization 
does not allow us to firmly establish whether it is one or the other. In the right-hand frame we show 
how this critical temperature depends on the magnetic field B. The dependence is monotonically, 
but non-linearly increasing. A similar behaviour was observed for h = 0. 

Our first main result is then shown in Fig. 5, displaying the phase diagram in the ji—T plane in 
the chiral limit, for a number of different value of the magnetic field B. The dashed lines show the 
part of the diagram where the transition is a cross-over, and the full lines where it is first order. 
At small chemical potential, increasing B corresponds to increasing T c . But interestingly, the lines 




Figure 4. Left: Minimum <j> of the effective potential for B = fj, = at the physical point as a function of 
temperature T. Ri£ 
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Figure 5. Phase diagram in the \i-T plane in the chiral limit for different values of the magnetic field. The 
arrows indicate the direction of increasing B. 

cross as the large- .B curves dip down at much smaller /j>b than the small- B curves. This means that 
at \x > 225 MeV increasing B at fixed fx corresponds to decreasing T c . 

At very large chemical potential, the phase transition becomes first order at a critical point (see 
inset). The crosses indicate the position of these critical points 5 (T*,/i*). Since these appear in the 
range [i > 225 MeV, the curves are ordered such that T* is an increasing function of B and fi* is a 
decreasing function of B. 

Additional features in the phase diagram have been reported in [31] at B = 0, h = and large 
chemical potential, using slightly different parameters. In this region the authors found that the 
first-order transition forks into two phase transitions. The left transition line is always of first order 

5 Note that T c is the temperature at a given /j,b where there is a transition, but is not the temperature at the critical 
end point. 
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Figure 6. Phase diagram in the ll—T plane at the physical point for different values of the magnetic field. 
The arrows indicate the direction of increasing B. 



while the chiral restoration transition line around the splitting point is initially of first order and 
turns into a second order transition for lower temperatures. 

Although we do see some possible additional structure appearing there, for instance in the \qB\ = 
0.2A 2 curve, this region seems to be quite sensitive to the discretization of the numerics, and we feel 
unable to make any firm statements. 

In Fig. 6, we show the phase diagram in the ll-T plane, but now at the physical point h > 0, 
again for different values of the magnetic field B. Again, the curves cross, now at a value of around 
ll = 210 MeV. In the inset, the crosses again indicate the position of the critical points. 

4 Discussion and conclusion 

In the present paper, we have mapped out the phase diagram of the quark-meson model in the 
presence of constant external magnetic field using the functional renormalization group. For small 
values of the chemical potential, our study shows that the critical temperature is an increasing 
function of the magnetic field B. However, for larger values of ll the curves for different B cross in 
the T — li plane and so the critical temperature becomes a decreasing function of chemical potential 
for large values of ll. This total inversion of the ordering of the curves takes place within a small 
range of chemical potentials of about 20 MeV, centered around 225 and 210 MeV for h = and 
h > 0, respectively. This behavior has also been seen in the NJL calculations of Refs. [14, 40] 
and the calculations of Ref. [41] based on holography. The phenomenon was appropriately dubbed 
"inverse magnetic catalysis" (IMC). Basically, there are two competing effects at work. Increasing 
the magnetic field B, effectively increases the particle-particle coupling and thereby enhances chiral 
symmetry breaking [34, 35]. The chemical potential ll tends to split the particles and antiparticles 
and thus works against chiral symmetry breaking. This is clearly demonstrated in the left half of 



Figs. 5 and 6; namely that for fixed n, the critical temperature increases with B, and for fixed B, 
the critical temperature decreases with \x. However, for larger values of B, the critical lines are 
crossing and so the critical temperature decreases with B (and still decreases with \i as well). Thus 
the effects of increasing the magnetic field must have changed qualitatively - namely that the cost 
of chiral symmetry breaking exceeds the gain [41]. 

Another main focus of our study was to pin down the critical endpoint, and its dependence on 
B. It turned out that the first-order line extends to higher temperatures for larger B; and because 
of the overall order- inversion of the transition curves the critical [i decreases with larger B. In 
particular, the largest magnetic field considered here corresponds to \qB\ ~ 5m 2 , a value expected 
to be achievable in heavy-ion collisions. In this case, the critical endpoint should be expected at 
(T*,(i*) ~ (81.4,215) MeV, rather than the B = values of (56.6,243). In particular, this means 
that the first-order region can be reached at about 10-15 percent smaller baryon density. 

In Ref. [21], the author have determined the phase diagram of QCD in the fx—T plane for B = 
using lattice simulations with physical quark masses and a Taylor-expansion technique around \i = 0. 
The prediction for the critical line is in reasonable agreement with the present results for small values 
of the chemical potential but the curves diverge at around /x = 150 MeV. Of course, for some (large) 
value of the baryon chemical potential, the lattice results can no longer be trusted so the discrepancy 
is not necessarily worrisome. 

D'Elia et al have carried out lattice simulations in a constant magnetic background at zero chemical 
potential [22]. They explored various constituent quark masses corresponding to a pion mass of 
200 — 480 MeV and different magnetic fields, up to \qB\ ~ 20 m 2 for the lightest quark masses. For 
these values of the pion mass, they found that there is a slight increase in the critical temperature T c 
for the chiral transition. These results have been confirmed by Bali et al [24, 25]. The same group 
has alo carried out lattice simulations for physical values of the pion mass, i.e. m n = 140 MeV. 
Their results which are extrapolated to the continuum limit show that the critical temperature is 
a decreasing function of the magnetic field [24, 25]. Hence the critical temperature for fixed \qB\ 
as a function of the quark mass is nontrivial. Moreover, at \qB\ = 1 (GeV) 2 the transition is still 
a crossover, although somewhat stronger than at B = This is in stark contrast to most model 
calculations that imply an increasing critical temperature as a functions of B. This is irrespective 
of whether one goes beyond mean field or not. Since this discrepancy is not understood, more work 
is needed to resolve the problem. 
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A Flow equation 



In this appendix, we briefly discuss the derivation of the flow equation (2.5). The starting point is 
the exact flow equation for the fc-dependent effective action Tk[(f>] [28] 



d k r k [4>] = ^Tr 



d k R k (q) T^+R k {p) 



(A.l) 



where the superscript n on T k [(f)] means the nth functional derivative of T k [(j)] and the trace is over 
the spacetime momenta p, and indices of the inverse propagator matrix. The function R k (p) is a 
regulator and is introduced in order to implement the renormalization group ideas: R k {p) is large 
for p < k and small for p > k whenever < k < A, and R\(p) = oo. These properties ensure that 
the modes below k are heavy and decouple, and only the modes between k and the UV cutoff A 
are light and integrated out. In order to proceed, we need to chose one regulator for the bosonic 
fields and one for the fermionic fields. The choice of regulator has been discussed extensively in the 
literature and in the present work we will use [42, 43] 

fl?(p) = {k 2 - p 2 )0(fc 2 - p 2 ) , (A.2) 

rFm - [i &VwX* - ') <* + M * - p2) ■ (A - 3) 

This regulator is very convenient in practical calculations since one can carry out the integral over 
three-momenta exactly. This turns the flow equation into a partial differential equation as we show 
below. 

The derivate expansion of the full quantum effective action reads 

m = [drf d*x {iz« [(V.) 2 + (Vvr) 2 ] + \zf [(doa f + (d 7r) 2 ] 

+... + ...U k (4>) + 4 3) ^[7o9o - 74/# + Z^inM + 9ki> [a - ij 5 r • vr] ^ + •••} , (A.4) 

where Zj?' are wavefunction renormalization constants, U k {4>) ls the scale-dependent effective poten- 
tial, and the ellipsis denote higher-order derivative terms that satisfy the symmetries of the effective 

(i) 

action. Note that the Z k are in principle different since Lorentz invariance is broken due to both 
finite temperature and finite density. In the local-potential approximation, one neglects the scale 
dependence of Z k and retain their tree-level values Z\ ' = 1. The flow equation then reduces to a 
coupled equation for the effective potential U k (<j>) and the Yukawa coupling g k . If we furthermore 
ignore the running of the latter, we obtain a tractable equation for the former. Ignoring the running 
of the Yukawa coupling, implies that the approximation does not correspond to the local-potential 
approximation [31]. 

The effective potential U k (<f)) is defined by evaluating the effective action T k [(p] for space-time 
independent fields and dividing by the VT, where V is the volume: 

u k {(f> un i) = -^r fc [0 uni ] . (A.5) 



In the following, we neglect the subscript uni. Inserting Eq. (A. 5) into Eq. (A.l), we obtain the flow 
equation for U. This yields 
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where <& = (a, 7r), i,j = 1,2,3,4, and the trace is over internal indices only. The symbols an d 
JJpj are defined by 
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Using the chain rule, the matrix appearing in the first term in Eq. (A.6) 
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where U' k = with p = \4> 2 ■ We notice that the matrix Eq. (A.9) is simply the inverse tree-level 
propagator if we make the substitution Vq — > Uk + R k (p). Integration over p gives 
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Summing over the Matsubara frequencies, we obtain 
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where = ^k 2 + U' k , u 2 ,k = ^Jk 2 + U' k + <f) 2 U k , uo q , k = y/tf + g^, N B {x) = \/{eP x - 1). and 
Np{x) = l/(e^ x± ^ + l). 

We next briefly consider the flow equation in a constant magnetic field B [10]. As we pointed out 
above, the matrices appearing on the right-hand side of the flow equation are essentially the inverse 
tree- level propagators. We can therefore find the flow equation in a constant magnetic background 



since the tree-level propagators are expressed in terms of the well-known solutions to the Klein- 
Gordon and Dirac equations in constant B field. The spectra are given by E 2 = p 2 + (2m + l)\qB\ 
for bosons and E 2 = p 2 + (2m + 1 — s)\qfB\ for fermions, where q z is the z-component of the three- 
momentum, m is the mth Landau level and s = ±1 is the spin variable. We therefore make the 
substitutions q 2 -» q 2 + (2m + and q 2 + (2m + l — s)\qfB\ in the regulators (A. 2) and (A. 3). In 

the case of the charged particles the sum-integral is replaced by a sum over Matsubara frequencies 
Pq = 2TrnT, a sum over Landau levels m, and an integral over momenta p z in d = 1 dimension: 



i 



\qB\T 



2ir 



E E 

P =2mrT, m=0 

\qB\T 
2ir 



dp z 
2vr 



E E 



dp z 
2tt 



(A.12) 
(A.13) 



P = (2n+l)irT, m=0 ' 

After integrating over p z and summing over Matsubara frequencies, one obtains. 

k 4 



dkUk 



12vr 2 



(1 + 2n B {u\ )k ) + (1 + 2n B (u 2t k) 



oo 



^ m=0 



2vr 2 ^ 



s,f,rn=0 



\qtB\k I 

\jk 2 - p\(qf, m, s) 9 (k 2 - p 2 ± (qf, m, s)) [l - ri£(u qf)k ) - , 

(A.14) 



where we have defined p\(q, m, s) = (2m + 1 — We close the Appendix by taking the limit 

B — > in Eq. (A.14). We change variable p\ = 2\qB\m, which yields p±dp± = \qB\dm. Replacing 
the sum by an integral, we obtain 



dkU k 



12vr 2 



(l + 2n B (w 1| fc) + 



<^2,fc 



[1 + 2n B (uJ2,k)) 



+ 



k f 00 I 

a / dp ± p ± Jk 2 -p 2 ± 9(k 2 -p 2 ± )[l + 2n B (io hk )} 

l,k^ Jo v 

J dp±P±^/k 2 -p\6 (k 2 - p\) [1 - n+ (w fe ) - np{u qjk )] . (A.15) 



2w 

N c N f k 

^g,k^ 2 Jo 



Finally, integrating over p±, the flow equation (A.15) reduces to Eq. (A. 11) 
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